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Abstract 

A spectrum of cosmic rays within energy range 10^^ — 3 x 10^^ eV was derived from the data of the 
small Cherenkov setup, which is a part of the Yakutsk complex EAS array. In this, work a new series of 
observation is covered. These observations lasted from 2000 till 2010 and resulted in increased number of 
registered events within interval 10^^ — 10^*^ eV, which in turn made it possible to reproduce cosmic ray 
spectrum in this energy domain with better precision. A sign of a thin structure is observed in the shape 
of the spectrum. It could be related to the escape of heavy nuclei from our Galaxy. Cosmic ray mass 
composition was obtained for the energy region 10^^ — 10^® eV. A joint analysis of spectrum and mass 
composition of cosmic rays was performed. Obtained results are considered in the context of theoretical 
computations that were performed with the use of hypothesis of galactic and meta-galactic origin of cosmic 
rays. 
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I. INTRODUCTION 



Energy spectrum of cosmic rays (CR) in energy range 3 x (10^^ — 10^^) eV could not be studied 
in detail with compact arrays due to their small acceptance at energy above 10^'' eV. At the same 
time this area of the spectrum is of a great interest, since local irregularities are manifested there: 
production of kinks (thin structure at 3 x 10^^ — 10^^ eV) arising from non-uniform distribution of 
heavier CR components in our Galaxy. On the other hand, this effect is smoothed by addition of 
a new component (of meta-galactic or other origin) to the cosmic ray flux near Earth. As a result, 
presence/absence of significant irregularities in spectra measured by various compact arrays allows 
one to speculate on the CR origin and propagation in our Galaxy [l,^. 

The Yakutsk array in this sense appears as a unique scientific tool. It is related to medium-sized 
arrays, capable of effective measuring of cosmic rays flux in a wide energy range (10^^ — 10^^ eV). 
Other important traits of the array are its model-independent technique of energy estimation of ex- 
tensive air showers (EAS) and the ability to track longitudinal EAS development by detecting the 
Cherenkov light emission. Factors mentioned above enable adopting the unique method, combin- 
ing the studies of CR spectrum and mass composition aimed at exploration of astrophysical aspect 
of cosmic rays |3, ^. 



II. METHODICAL ISSUES 

For more than 15 years the small Cherenkov setup has been operating as a part of the Yakutsk 
array. It measures Cherenkov light emission from EAS of lower energies (see Fig{T]) using standard 
detectors which are designed to operate in winter conditions. The area of modern prototype was 
significantly increased in comparison with the original setup, its border forms a circle of 500 m 
radius. The number of optical detectors was also increased (see FigH]). Table H] presents the 
information on operation of the setup (on annual basis) combined with mean spectral atmosphere 
transparency at wavelength 430 nm. 

All the information on each shower is stored in the database, which is controlled by a software 
complex. This program includes units for gathering, sorting and storing of the experimental data. 
It also includes mathematical units for data processing and statistical analysis. The results of the 
analysis are presented below. 
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FIG. 1. Location of observational points of the small Cherenkov setup. 

TABLE I. Characteristics of different observational periods 
Years 95 95-96 96-97 97-98 98-99 99-00 00-01 01-02 



T, min 9960 19164 22563 26586 28457 29683 26530 26947 
iV, events 10956 18973 27076 26055 38418 27161 28652 25761 

rel. units 0.61 0.59 0.62 0.58 0.64 0.60 0.62 0.61 
Years 02-03 03-04 04-05 05-06 06-07 07-08 08-09 09-10 
T, min 26912 33922 29613 32044 32640 29725 25514 35609 
iV, events 23548 31718 30649 30634 34272 31211 25131 44867 
p, rel. units 0.59 0.61 0.63 0.61 0.63 0.63 0.61 0.65 



A. Selection, processing and choice of the classification parameter for showers 

In order to reconstruct cosmic ray spectrum we used following criteria for shower selection: 
1. atmospheric transmittance px > 0.65; 2. shower axis must lie within 250 m from the center 
of array (for showers with primary energy Eq < 3 x 10^^ eV) and within 500 m (for showers 
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with Eq > 10^^ eV); 3. zenith angle of shower arrival direction 9 < 50°; 4. chance probability of 
shower detection Wr > 0.9. Introduction of these criteria was mainly defined by the parameters 
of detectors used in the experiment (mostly by aperture and thresholds of Cherenkov light detec- 
tors) and by atmospheric conditions during optical observation. Since shower events were selected 
by different triggers, during transition from one trigger to another a threshold effect appears. It 
manifests itself by sudden local increase of intensity in the spectrum. According to our calcula- 
tions the magnitude of this effect can reach ~ 30 %. Technically, the softening of this effect has 
been achieved by extension of the area controlled by the trigger-50 and by its overlapping with 
other triggers. For example when triggers coincided, having been activated simultaneously by the 
same shower. In this case, the number of skipped showers was minimized. Another method is 
introduction of corrections obtained in simulation of measurement. 

According to the criteria described above, the data bank of showers was formed by the param- 
eter (5(150), a density of Cherenkov light flux at r = 150 m from shower axis. This parameter 
was derived from readings of Cherenkov detectors located within 80 — 200 m from shower axis. 
The structure of the data bank was defined by the task — it was formed strictly by those periods of 
observation, which confirmed to adopted shower selection criteria mentioned above. We suppose 
that chosen conditions are sufficient to avoid distortions related to experimental and methodical 
errors in reconstruction of cosmic rays energy spectrum. 



B. Monitoring of the atmosphere 

It is believed that photon losses in clear atmosphere arise from Relay scattering (5 % from total 
flux). In real conditions there is significant loss in received light due to Mi-aerosol of various size. 
In winter (in the region where array is located the climate is sharp-continental) the atmosphere 
above the array is non-standard, its parameters change significantly from autumn to winter and 
vice-versa. According to work [|5D all this factors should be tracked on an operational basis and 
taken into consideration when analysing different observational periods. On Fig. [2] perennial data 
are presented on Cherenkov light transmission in atmosphere during different periods of optical 
observations. These data were utilized during generation of shower samples from which cosmic 
ray spectrum had been calculated. 
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FIG. 2. Atmospheric transmittance averaged by months for wavelength 430 nm. 



C. Estimation of primary energy 



Energy was estimated with quasi-calorimetric method, from joint measurements of main 
shower component. Recently, thanks to significant increase in acceptance for showers with energy 
< 500 PeV, it has become possible for Yakutsk array to measure total number of charged particles 
on observation level with the precision (20 -ir 30) % and density of Cherenkov light flux at any 
given distance from shower axis with the precision 15%. Such distances for small and large 
Cherenkov setups are r = 150 m and r = 400 m accordingly. The densities of Cherenkov light 
fluxes at these distances were adopted as classification parameters. On fig. |3l energy dependence 
of classification parameters is shown. From the data presented on Fig. [3l a connection between 
classification parameters and primary energy of a shower Eq has been derived: 



Eo = (9.12 ±2.28) X 10 
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FIG. 3. Energy dependence of classification parameters. 

Expression ([T]) is used in energy range (5 -r- 500) PeV and expression ^ — to estimate primary 
energy of showers above 5 x 10^^ eV. 

D. Calculation of the Seff. 



The effective area (Sgff.) of the array was taken into consideration for each energy interval ac- 
cording to the method of shower selection. If time of observations and the spatial angle are always 
defined accurately, then underestimation and overestimation of the Sen. usually results in system- 
atic errors in the estimation of the absolute EAS spectrum intensity. This is why the thresholds of 
Cherenkov detectors (3 ■ 10^ photons/m^) and the effectiveness of shower selection system were 
considered during the estimation of the Sgff., i.e. various configurations of triggers of the small 
setup. For the small Cherenkov array such trigger is — a triple coincidence from Cherenkov 
detectors. In the latter case a triple excess of the signal amplitude over the night sky background is 
required. The calculation had been carried out with the use of Monte Carlo method within energy 
Eq intervals with logarithmic step logio = 15.1. For each step 1000 artificial showers were 
thrown with axes distributed within a circle of 250 m radius for showers of a smaller energies and 




FIG. 4. Effective area Sgff. of the small Cherenkov setup for showers from different triggers. Solid line — 
showers with smaller energies (axes within r = 250 m), dotted line — showers with Eq > 10^^ eV (axes 
within r = 500 m). 

500 m for showers with > 10^^ eV (see Fig. H]). A transition from classification parameters 
to the energy was made according to empirical formula ([T]). When calculating intensity of the 
spectrum, showers were considered with 9 < 50°. The time of the observation is presented in the 
Table E 



III. RESULTS 

A. Energy spectrum of cosmic rays within energy interval 10^^ — 10^*^ eV 

Energy spectrum of cosmic rays with the account of the new data on showers from the selection 
is presented on Fig. [51161 Statistical accuracy allows one to speculate of a thin structure within the 
shape of the spectrum. A comparison with model calculations performed within various hypothe- 
ses of CR sources and models of their propagation in the galaxy gives possible interpretation of 
the experimental data. 
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FIG. 5. Interpretation of the CR spectrum according to scenario 1. The Yakutsk data are denoted with red 
triangles, black squares and blue stars denote results obtained at KASKADE Grande. Solid line — part of 
the spectrum formed by galactic CR, dotted line — extragalactic component, dashed line — a component 
yet to be explained. 



Here we consider two possible scenarios of generation of the spectrum. In the scenario 1, the 
kink at E = 3^^ eV and subsequent increase result from galactic component (up to 10^^ eV). In 
the region (~ 3 x 10^'' — 3^^) eV, the spectrum is shaped by unknown component; here, one may 
speculate of CR particles interaction with galactic wind and shock acceleration. In the scenario 2, 
the galactic component extends to 3^^ eV thanks to acceleration by supernova remnants M>[B\- This 
model shows good agreement with experimental data and helps to interpret observed spectrum in 
the range 10^^ - 3 x 10^*^ eV. 

As it is seen from figures, within the energy range (5 — 8) x 10^^ eV there is a small peak 
generated by iron nuclei. At lesser energies, where a peak from the CNO group is expected 
according to the model, a slight increase in the intensity is observed. According to our data and 
the data from KASCADE-Grande, this exceeding is not significant and could be related to the 
presence of CR of another origin in the total flux [|20. 
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FIG. 6. Interpretation of the CR spectrum according to scenario 2. Designation of the experimental data 
is the same as the previous figure. SoUd curve — theoretical predictions obtained by Berezhko and Volk 
(2009), dotted curve — by Ptuskin and Zirakashvih (2010). 

B. Mass composition 

The mean natural logarithm of the CR atomic number (In A) derived from the Xmax measured at 
the Yakutsk array is shown on the Fig. |71 The (In A) calculations were performed with the use of 
the Xmax predictions of QGS Jet II and SIBYLL models for proton and iron and with the relation dS]) 
proposed by Horandel et al 

(lnA) = ^^^i:i4^.1n56 (3) 

•^max -^max 

The data from the Fig. |7] point towards the slight change of the (In A) value right after the 
kink in the spectrum. For example within the energy range 10^^ — 10^ eV the (In A) value gets its 
maximum ~ 3 at i^o = (5 — 8) x 10^^ eV (the composition becomes heavier) and above the energy 
3 X 10^ eV a decrease of the (In A) is pointed out (i.e. the composition becomes lighter). The data 
from other FAS arrays testify of the same tendency [8]. If one plots the mass composition data 
versus the CR energy spectrum (see Fig. IT]), then the coincidence between the peak of CR intensity 
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FIG. 7. The dependence of the (In A) from EAS energy. Predictions of the QGSJet 11-03 model are 
presented with circles, squares represent mass composition estimation according to the SIBYLL2.1 model. 
Curve — expected composition from supernova remnants (Berezhko and Volk, 2009). 

and the maximum of the (In A) value could be clearly seen. Then it follows that the nature of the 
peak in the spectrum is related to a heavier component of cosmic rays. 
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